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An approximate analytical queuing model for the intraship transfer 
of material phase of the replenishment at sea problem is formulated. 
The physical model has two parallel sets of three service facilities 
in series that move material to a final stage of three service 
facilities in parallel. Exponential service times and finite queue 
Space are assumed for each facility, except that the initial facilities 
always have infinite queues. This implies that the output from an 
individual facility is approximately Poisson. The mean output rate for 
the queuing system is obtained. Several parametric studies on critical 
parameters are performed, with the results presented in graphical format. 
The model is easily enlarged to allow a general number of facilities in 
series in any of the parallel paths. Generalization to more than two 


paths in parallel for the initial stages may be possible. 
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I. INTRODUCTION 


Underway replenishment (UNREP) of the ships of the United States 
Navy at sea is a mandatory requirement in maintaining a continuous 
posture of readiness to meet the present and future global commitments 
and responsibilities of this country. However, replenishment of ships 
at sea should detract as little as possible from the primary missions 
of a fleet. This implies that the goal of the UNREP operation is 
complete resupply of combatant ships in the minimum possible time. 

There are several phases of the replenishment operation, | possibly 
not all distinct. One list of phases might be: Long range planning 
required to schedule a replenishment both in terms of time and place; 
the actual movement of ships toward the prescheduled rendezvous; when 
several supply and combatant ships are involved in a replenishment 
operation, the scheduling and movement of all the combatant ships 
through the replenishment cycle; and the actual transfer of material 
from the hold of the supply ship to designated storerooms on the 
combatant. It is obvious that all phases require careful attention 
if the operation is to be efficient. Consequently, all phases are 
fruitful study areas in a continuing effort to improve UNREP efficiency. 

The study reported herein dealt with the actual transfer operation. 
Perhaps it is the most important phase in that the ships involved are 


most vulnerable, their capability to maneuver reduced and their 


| t Throughout the paper the terms UNREP and replenishment operation 
will imply the entire operation; planning, scheduling, shipmovement, 
etc., while the term transfer operation denotes the movement of material 
from hold of the supply ship to storeroom of the combatant. 


watertight boundaries open to allow passage of material. Unfortunately 
in the past it appears that this phase has received less attention than 
it warrants both from the standpoint of detailed study and the introduc- 
tion of improved material handling equipment and techniques. 

Recent efforts have been made that should considerably improve the 
overall speed and efficiency of the transfer operation. Improvements 
include the introduction of new supply ships with the capability to 
employ the unit load system of material storage and transfer, more and 
better mechanized material handling equipment on both supply and 
combatant ships, and a realization of the problem by ship designers 
which has resulted in improved match of replenishment stations between 
supply and combatant ships. The unit load concept provides for mate- 
rial to be stored on the supply ship in the configuration for future 
transfer. Full implementation of the concept will maintain the unit 
load intact through ultimate storage in the storeroom or magazine of 
the combatant receiver. Although this concept may never be completely 
feasible, especially in the case of smaller combatant ships, it should 
provide a vast improvement to perhaps the most important and demanding 
of all UNREP operations, the rearming of aircraft carriers. 

The transfer operation can be represented as a queuing problem. 

It is a system of parallel queues in series with movement between the 
parallel paths. Attempts at reaching an analytic solution are further 
thwarted by the finite waiting room in the queues involved. Formal 
efforts to analyze the problem have been few in number. The powerful 
computer simulation tool has been the general mode of analysis. An 
early simulation was conducted by Stanford Research Institute [8]. A 


detailed computer simulation study has also been done by Alfred M., 


Feiler of Project TRANSIM, University of California at Los Angeles, 
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Department of Engineering. Computer simulation studies are currently 
being done at Stanford Research Institute and Center for Naval 


Analyses, 


It, BACKGROUND 


A. THE TRANSFER OPERATION 

Although some transfer of material between ships at sea is 
presently accomplished by helicopter, the vast majority of material 
moves between ships via attached horizontal rigs. A typical move- 
ment sequence for material transferred at sea would commence with a 
vertical lift from the hold of the supply ship to the main or transfer 
deck by means of an elevator or a boom and windlass. Once on the 
transfer deck the material will be moved to the transfer station by 
fork lift truck or hand pallet truck. The intraship transfer of 
material is accomplished by one of several different cable and pulley 
arrangements. The type employed on a given transfer will depend upon 
the types of ships involved and the amount and type of material to be 
transferred. Material landed on the transfer station of the combatant 
ship is cleared from the area by fork lift or pallet truck, Eventually 
the transferred material will be transported vertically down to its 
designated storeroom or magazine by elevator. 

The type of movement just discussed employed the unit load concept 
required by the model developed herein. As noted earlier, many small 
combatant ships must break up the unit loads into smaller units for 
movement beyond the transfer station. Many times hand movement is 
required because of space limitations in the ship's internal passages. 


The approximate model developed by this study could still be used for 


this type of movement using a mean service time for the unit load 
transferred by pieces. 

Queues of material and subsequent blockages of material can and 
do form at several points along this path. On the supply ship, a 
queue of material will probably exist on the transfer deck near the 
point where vertical movement from the hold is completed. A queue 
of unit loads will also exist at the transfer station as material 
awaits movement to the combatant. The transfer station on the 
receiving ship is a critical blocking point as transfer of a load 
from the supply ship to the combatant may not begin until the landing 
area has been cleared of the previous load, Generally, material 
cleared from the landing area on the combatant ship is temporarily 
stored in what is referred to as a marshalling area. Material in 
the marshalling area is distributed to a strikedown route and then 
to eventual storage in the assigned storeroom., In many cases, 
additional queues of material will form along this strikedown route, 
especially at the top and bottom of vertical movement leg(s). 

Much of the material transferred at sea presently is stored as an 
individual component or packaged with similar items in such a manner 
that hand movement is possible. These individual components are 
combined into loads that can be conveniently handled by the transfer 
rig between the two ships. This load makeup is accomplished either 
in the hold or on the main deck of the supply ship. The process is 
essentially reversed on the receiving ship with the breakup of loads 
commenced early in the storage routine, especially on smaller ships, 
where it is generally cleared from the landing area piece by piece. 


Introduction of more and better handling equipment on all ships and 
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movement of material via unit loads through all or most of the cycle 
will vastly improve the efficiency of the operation. This step forward, 
although far from universal reality, is in progress, 

Careful consideration indicates that the three distinct entities 
having a major effect on the movement of material through the system 
are the service or handling rates through all phases, the storage 
capacity of the various queues that inherently form and the origination 
and final distribution of the transferred material. The importance of 
the handling rates is obvious to the extent that discussion here is un- 
necessary. However, the importance of the queuing (temporary storage) 
space and the distribution of material requires some explanation. Any 
amount of free open space on any ship is always scarce, this includes 
space for temporary storage of transferred material. Perhaps the single 
exception is the hanger deck on an aircraft carrier, (aircraft carrier 
transfer stations are located off the hanger deck), however, even here 
Space may be at a premium as it is normally taken up by aircraft in 
various states of overhaul. Therefore, the transfer process is often 
slowed because of blockages caused by the inherent space limitations 
on board ship. 

Limited space on ships is also partially responsible for transfer 
slowdowns that can be attributed to excessive distribution to a single 
material handling facility, e.g., a vertical strikedown route on the 
combatant ship. Space available on a ship restricts both the size of 
storerooms and vertical and horizontal passages leading to the store- 
rooms. In order to provide for convenient and rapid retrieval of a 
given item, it must be stored in a predesignated place. Further, 


many items require special storage facilities that are available in 


11 


only a single location. These special facilities can range from 
refrigerated storage for fresh and frozen provisions to storage trays 
and racks for most ordnance. Even if the above difficulties could be 
overcome, space limitations again would be limiting in that storerooms 
are generally kept full so that the only storage space available is 
the designated one. Therefore, it is often impossible to change 
material destination to accommodate transfer efficiency. A possible 
solution to blockage on a route would be an alternate route for some 


of the material. Unfortunately, alternate routes may not always exist. 


B. PARALLEL/SERIES QUEUING MODELS 

Solutions for a single stage of several parallel queues, with 
various types of input and service times, have been used for some time. 
Knowledge of single server queues in series is more limited and is 
generally restricted to Poisson input to the first stage and exponen- 
tial service times in all stages with infinite queues allowed. 
Information in the literature of a combination of queues in series 
and parallel is virtually nonexistent. 

Saaty [ 6| gives a survey in the area of series queues. The 
bibliography for this work is nearly complete and he has recently 
addended the bibliography to include recent contributions in the 
field ad Burke Ly reported the fundamental result for queues in 
series with Poisson input to the first stage, exponential service 
time in all stages with infinite queues allowed, showing that the 
output from each stage and therefore the input to the next stage 
is also Poisson, Hunt [4] using the assumption of exponential service 
times did initial work in the area of allowing only finite queues to 


form before each series facility. Hillier and Boling [2] extended 


b2 


Hunt's work under the exponential service time assumption and also 
solved the problem for Erlang service times for each stage. 

Hillier and Boling concurrently reported an approximate procedure 
for determining the mean output rate for a system of exponential 
service time facilities in series with finite queues allowed before 
each facility except the initial stage where an infinite queue was 
assumed to exist. Hillier and Boling compared results obtained with 
the approximate procedure with the limited exact results obtained. 
The comparison showed very little error in all cases with the magnitude 
of the error decreasing with increased allowable queue size. 

The Hillier and Boling approximate procedure has the distinct 
advantage of being computationally feasible, even for very large 
problems, This approximate procedure formed the basis for the study 
reported here and it will be discussed and outlined in some detail in 


the next section. 


C. APPROXIMATE PROCEDURE 

The important features of the Hillier and Boling [ 2] procedure, 
for approximating efficiently the mean output rate (R) for a series 
of single server exponential service time facilities, may be outlined 
as discussed below. Assume that there are N facilities in series each 


and waiting 


3 3 


with an exponential service rate of A, ga WS Ae N 
imeomen. , 1 = 253, °°" 5 Nj where the M;'s include the item being served 
and that an infinite queue always exists before the initial facility. 
The basis for the procedure comes from the results due to Burke [1] 
which indicated that when the M,'s are equal to infinity then each 


service facility would have a Poisson input. Reich [5] showed further 


that in this case the queue sizes would be independent. This implies 
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that each facility could be analyzed separately using the standard 
M/M/1 model. The approximate procedure with finite queues allowed 

is then based on the assumption that the output from each stage is 
approximately Poisson. Each facility is analyzed individually by 
using the single-server queuing model M/M/1 with finite queue allowed. 
Note that R is merely the mean output rate for each of the facilities. 

The procedure then assumes that there exists an effective mean 
arrival rate AD (when there isn't blocking in the cape facrrimeyy 
to the ah facility. It is also assumed that the queue lengths before 
each facility are independent so that the unconditional probability 
(1-PG-1)) that the equ rye facility is busy is equal to the conditional 
probability that it is busy, given that the queue for the jth facility 
is not full. Thus it follows that, 

(1) AY = Hj) , if n&M, 
0 , n=M;+1 
where n is the number of customers in the queue for the jth facility 
(including any customer being held at the (jeaye" facility because of 
blocking) and j = 2,3, nm. 

The effective service rate for the jth facility, AC}t, will be 
somewhat less than £,, j=2,3, ‘++ ,N-l, because time spent in holding 
a customer after service due to blockage in the (j+1)S© queue must be 
included. Now the mean output rate of a single-server queuing system 


is the product of the mean service rate and the probability that the 


server is busy. Therefore, it follows that, 


Gis = _p(j) — soe 
(22) ACG) = r/a-eG)) —, for j = 2,3, °*° GN. 
However, for the final server the customer is always released immedi- 


ately, tvew, WD) - Ay, Therefore, 
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(3) R = Ace), 


Thus R may be determined if ve is known. 


For a single server queuing system with finite queue, it can be 
. [ 
shown that 3 | 


i+] 

| Cl) 
(4) es) eee 

2 M.+2 

l-y, J 


where for the case of the approximate system, 


= Fp Gp 
yy > Aes AA tt 
It is also clear that since the queue before the first facility is 


infinite, pir) = 0 or ap) = 1, 


It is now possible to determine the true value of R by first 
assuming a trial value of R, call it R*. Since -P is known, Fwd 
can be determined from equation (1). It then can be used along with R* 
in equations (2) and (4) to determine (1-P62)) Then 1-P(3)) , (1-P64)), 

; c-P?) can be obtained in a similar manner so that R can be 
calculated from (3). This calculated value of R is compared to R* and 
a new trial value assumed and the process repeated until the true value 
of R is obtained when R* - R = QO. 

Hillier and Boling have simplified the above procedure for solution 
by numerical methods and the equations used in this study are presented 


below: 


Combining equations (1) and (2) yields, 


| 
No 
Lo 
Z 

' 
i_— 


Zee, _ 


wel 
ij l- 0s [My 


ba i 


L5 


then, 


(5) - Pt)? = R Yas j ms Zi 5 avax vlad 


so that y can now be found by combining (4) and (5), as the positive 


real root of the equation, 


M1 
(6) (l-y.J) 
yj —— oe R - = 0 
M +2 (j=) 
-y (142 
j j-4 O 


Hillier and Boling show that this root exists and is unique for values 
* tt 
of R’ that are sufficiently close to R, the true mean output rate. 


Let, 


fou = a j ’ j a ee oe eel 
J oj 


and 
Be. a - e.{ ) == ais j= oe. 3., gaa , N-1 
See j° 7 “ee 


where R is a trial value, 


c, w__ 
i *”* 
AA; | j- 1 


’ 


* os 
and ~~ = £509), Y corresponding to the unique root of g(R,y,) = 0. 


4 
Finally as yields the calculated value of R where 
I 
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Tit, FORMULATION OF THE PROBLEM 


For the replenishment at sea transfer problem, as with many other 
physical problems that can be formulated in a queuing~ theory context, 
one of the interesting questions that can be asked is, "How long will 
the operation take?" This question was answered in the formulation 
for the UNREP transfer problem by determining the mean output rate (R) 
for the system. 

As previously noted, each at sea transfer of material will be 
somewhat a a queuing standpoint because of the ships 
involved, volume and type of material transferred, material origin 
and destination, types of material handling devices employed, etc. 
Thus, the number and size of queues formed during any transfer may 
differ considerably. 

A completely general model to represent the transfer at sea opera- 
tion would require up to five parallel paths from the supply ship to 
the combatant ship. On each of the parallel paths there are several 
material handling facilities with associated queues and waiting space 
on both the sending and receiving ships. Additionally, the general 
model would require movement of material from one parallel path to 
another at least once on each ship. 

The general model described above may be simplified considerably 
and still represent accurately many transfer operations. The model 
developed for this study utilized two parallel paths from the hold 


of the supply ship to a marshalling area on the combatant ship. From 


each of the two marshalling areas the material is moved to any one of 


three strikedown routes. Figure 1 is a pictorial representation of 


1b, 


this model. The inclusion of only two parallel paths was done for 
computational efficiency, noting also that the vast majority of at sea 
transfers are conducted with but two attached rigs. A generalization 
of the model for more than two paths is discussed in section VI. 
Restriction of movement between parallel paths on the supply ship is 
the other major simplification made for this study. It was felt that 
such an assumption was not overly restrictive and accurately represents 


many actual transfers. 
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Figure 1 


If individual material handling facilities in parallel are referred 
to as a stage, then the physical location of the first stage queue of 
the model can represent the hold of the supply ship. The second stage 
queues are formed by material on the main deck of the supply ship 
awaiting intraship transfer. The marshalling areas, which will gener- 
ally be adjacent to the transfer stations on the receiving ship, form 
the waiting space for the third stage queues. The fourth stage queues 


may be located physically at any one of several places on the individual 
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strikedown routes to the designated storerooms, e.g., on the transfer 
deck near an elevator or hatch, or at the bottom of a vertical ~ovement 
leg. 

As stated previously, service of material in the two third stage 
queues consists of movement to any one of the three strikedown routes, 
In the model, this movement is governed by two probability distribu- 
Beens™p, “arid q;j) i = 1,293, feorwthe upper and Wower routes respectmvely. 
The distributions are determined as that percent of the material moving 
across one of the initial routes that is scheduled for eventual storage 


th otriked@WnUroubela SlaemccueW auc see meen 


in a storeroom on the i 
effect this determination are the amount and type of material trans- 
ferred, the storage of material on the supply ship, the location of 
the storage space available on the combatant ship, and the availability 
of alternate routes to final storage. It is felt that the distribution 
approach as a means of representing the actual movement during transfer 
is reasonably flexible to cover accurately many transfer operations. A 
primary consideration in this regard is the fact that material is 
selected for movement in the most convenient manner, i.e., a load is 
selected for movement because it is available to the transfer facility 
or because it can be moved ahead of other unit loads to a nonblocked 
facility upstream, 

Within the model it is assumed that all material is moved as a 
unit load, Solution for smaller units would overly complicate the 
problem, however, the unit load assumption is enhanced by the present 
trend toward actual movement of material in this manner. This assump- 


tion had a further and probably a greater impact on obtaining an 


eventual solution because of the effect it had on the transfer rates, 


Ni 


Since all material moves in thisS manner it is reasonable to assume that 
a given facility can handle a unit load of 500 1b. bombs or 1000 1b. 
bombs or even fresh provisions, for example, at approximately the same 
rate. Thus the unit load concept helps justify the requirement that 
each facility must move all types of material at a fixed rate during 
each transfer investigation. There is of course some allowance for 
variation, e.g., if two types of material had different service rates 
for a given stage and each moved on different parallel paths, then the 
service rates for the two paths could be adjusted to represent more 
accurately the material being handled. 

Notation within the model is such that the first subscript (i) 


represents the gth 


h 


stage in series and the second subscript (j) repre- 
i : ; 

sents the j parallel element in a stage, e.g., Mh, is the service 

rate for the second parallel facility in the first stage and My] is the 


number of unit loads that can queue before the first parallel facility 


in the third stage. 


IV. PARALLEL SERIES SOLUTION 


The flow of material through the physical model of figure 1 may be 


represented as in figure 2. If R is the mean output rate, it follows 


Ry 
_ 
a. R 
_ 
SS RELI 


Figure 2 


20 


from previous discussion of the p,. and q; distributions and 
figure 2° that, 


Rares 
Z 


R = Ry =F Roy aF Riqt 


ee 


2... SpA ane 
tr Poy” 49% 


Tir Pa] ae 


where Ry > Ro » Ro; Rit and Rott are the mean output rates of their 
respective routes. For any UNREP that can be represented by this model, 
the value of the mean output rates will be determined by the various 
Service rates, queue waiting space, and the p; and q;, distributions, 
If the system were symetric, that is if the service rates and queue 
lengths for the facilities in each stage were the same and the Ps and 
qs distributions moved material symetrically to the final stage, then 
R, would equal Ro. A solution to this problem could be easily reached 
by the Hillier and Boling procedure by considering only half the system 
and solving for R/2. However, these conditions will not often hold and 
in general, AS will not equal Ro: 
When R and Ro are not equal, the Hillier and Boling approximate 
procedure cannot be applied directly to the system to obtain a solution 
for R. The problem is that not only must a trial value for R be 
assumed, but also values for Ry and R, must be assumed such that 
RI + R> = R’.. where the R“'s are trial values. Now there are infi- 
nitely many combinations of R = R* + R* = Re The problem thus 


i! 2 
reduces to finding the values of Rj and Ry that will yield the true 
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mean output rate. The true values of the R's will be complicated 
functions of ., M.., p., and q_. It will not, however, be necessarv 
ij ig i it 

to determine this function, as the solution may be reached by simply 
maximizing the sum R, + Ry for a given value of e 

Before proceeding to show how the true values of Ry, Ry and R may 
be obtained, the Hillier and Boling equations should be adapted to the 
problem, Assume that for a given problem the true values of Ry> Ro» 
and R are known. The solution for the output across the upper route 
may be found by applying the approximate procedure previously outlined. 


“ 
First Yo, may be found as the unique positive real root of 


* 

: x, = f Poe c a6 

86! l Yo? 91 791 21°21 
* * xe 
where Cc = R / o% a. = 1). Then, 
24 / ke] ( | ) 
*M 

: : ¥ (ih ory 
Ya ee 21 
“a 2 N 

21 a ‘= 20 \-y var" 


* x 
Ya1 and then» X.,, may be found in a similar manner. Also, X*, may be 
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found using the same technique using RA in place of Ro. Then since 


the fourth stage is the final stage “f" i= 1,2,3, can be determined 
directly as 
,* -* 
x et ( Ay X41? + MX 59%iy eGaia3 
i i hh = 


Mi 


* * x 
so that finally R = XA v KG,» + X od 


The solution procedure for the interesting case of unequal parallel 
de se 
outputs can now be outlined. As before, Ris assumed and them Ry and 
* # * * 
Ry are chosen such that R} + Ro = R . The problem is then solved as 
just outlined to obtain a computed value of the mean output rate, R. 


° * ¥ j ° a 
Using the same value of yA R; and R» are varied (maintaining the 
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relationship RY AP R> = R”) until the maximum computed value of R is 
found.) As in the original Hillier and Boling procedure, the computed 
value of R is compared to R™ and if R = ne the true mean output rate 
has been found. If R is not equal to R™, then a new value of R”™ is 
assumed and the procedure repeated until the true mean output rate is 
obtained. This procedure converges rapidly if two techniques are 
employed. When selecting a new value for Ras select R™ (new)=(R FR) /2. 
Also when commencing each iteration Ri (cal Boescia) = RY/R5. 
Justification for this procedure follows from physical flow laws. 
It is easy to see that the procedure is correct at the boundaries. For 
example, assume that the system was such that one of the 4; ;'s = O in 
the lower path. Then maximum R will be obtained when R® = RT and R5 = 0. 
A computer program to solve the problem of figure 1 is presented on 
page 34. The program was employed to obtain the results presented in 
section V. The wide range of input data was handled with no difficulty. 


Average computation time for the IBM 360 was 2 seconds. 


V. PRESENTATION OF RESULTS 


ee ae ee 


As noted previously, the three types of parameters within the 
approximate transfer model that have an effect on the mean output rate 
are the service rates of the material handling facilities (M4; 5's), the 
available waiting room in each queue (Mj;'s), and the distribution of 
material between the series facilities as governed by the pj and q,; 
distributions. The model can be employed to study UNREP transfer 
Operations in at least two ways. First of all, general parametric 
Studies can be made to determine the effect of varying the parameters 
in all possible cowhinations. Also, more specific studies may be made 


employing known handling rates, queuc space and expected material 


A 


distribution for specific ship types, to determine critical areas that 
can be improved to increase overall transfer efficiency. A second use 
for the model would be as an aid to the commander, perhaps even on as 
low a level as commanding officer of a supply ship, to help him in the 
planning for replenishment operations. All that would be required is 

a knowledge of material handling characteristics of the ships involved, 
the amount of material to be transferred and a small computer which is 
generally carried on board many supply ships today. 

The primary advantages of the approximate procedure over known 
computer simulation models are the simplicity of the analytical model, 
the inhePerd computing efficiency and the ease with which inputs into 
a computer program may be changed. Possible objections to the model 
Over computer simulation models would include the limitation that 
different types of material must be transferred through a facility at 
the same rate and that unit loads of each type of material would occupy 
the same volume or area of queue space. The exponential service time 
assumption with its memoryless property, i.e., the probability of a 
unit completing service in a time interval At is constant, no matter 
how much time has elapsed since the last completion of service, also 
may partially objectionable. A final objection may be raised over 
the final output of the model which is a Steady state result, i.e., it 
is assumed that the system has been in operation for a long time and 
is in equilibrium. This objection is at least partially overcome by 
the generally large volume of material transferred during the UNREP's 
that are of interest and the comparatively small queue sizes allowed, 
thus assuring that equilibrium conditions will be nearly attained 


during an operation. 
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Actual UNREP transfer data was not available in the format required 
for entry in the model. Therefore, numerical results were necessarily 
confined to several parametric studies. Graphs of the study results, 
along with the input data for the constant parameters, is presented in 
Appendix A. 

Parametric studies 1, 2 and 3 all vary the parameter Py> distribu- 
tion of material from route Ry to route R,, over the range 0.333 to 1.0. 
The parametric variation in the three studies involves the qd. distribu- 


Eton as follows: 


Study one (page 29) - - - dee = Wes. i132) 12 a6 
Study two (page 30) - - - q, = OZ. 1, = Lao 4 
Study three (page 31) - - Fadl 0. 4e age = q,=0 3 


All three studies indicate the marked decrease in overall transfer 
rate as congestion at the upper strikedown route (facility (4,1)) 
increases, As predictable, the decrease is most rapid for study 3 
where the flow to facility (4,1) and, therefore, the congestion at the 
facility is greatest. An additional effect, evident in all three studies, 
is the difference between Ry and as that can be attributed to the slight 
reduction in handling rates of facilities (2,2) and (3,2) and the 
smaller queue allowed at facility (3,2). 

Parametric study 4 gives a good indication of the efficiency that 
can be gained by increasing the handling rates of the facilities in a 
single stage. In this case, the service rates of facilities (2,1) and 
(2,2) were varied together over the range 0.5 to 1.0 loads/min, with 
all other service rates fixed at 0.5 loads/min. The graph of the study 
(page 32) indicates that small increases above the system average will 


have a significant effect. Study 5 varied the queue space in facilities 


7s, 


(3,1) and (3,2) over the range 5 to 10 loads with the remaining queue 
spaces held fixed at 5 unit loads. Once again a graph of the results 
(page 33) indicates that small changes above the norm are worthwhile. 
Comparison of studies 4 and 5 shows that increased service rates have 
a greater effect on efficiency than increases in allowable queue size. 
This is understandable and also encouraging, as generally, space on 
ships will not be available to increase temporary storage to promote 


increased transfer efficiency. 


VI. BAITENSTON OF THE. MODEL 


There are several ways that the specific model developed here can 
be extended. Unfortunately, the most interesting extension becomes 
computationally difficult. 

The most obvious extension follows directly from the Hillier and 
Boling {2 | series development where the number of stages in series is 
N. The parallel series model developed here may be similarly extended. 
Referring to figure 2, this would mean that there could be N and ~ 
separate service facilities in the two initial parallel routes and Ny» 
N and N Separate service facilities in the three final routes, 


II’ IIT 


where N N N and N are only restricted to be finite. 


al iS | ali lus III 
Solution of this generalization should present no special difficulty. 
The interesting extension fron an UNREP point of view would seem to 
be that which generalizes the number of parallel paths allowed. Such a 
generalization is represented in figure 3. 
In solving such a system, it is first of all evident that if R > 


R R_ are known, then R R R,, are uniquely determined 


ur, r=" 


by the Pi, distributions. Where Pa j is the distribution of material 


26 


from the aa initial stage to the men final stage, Since it must hold 


that 


and 


JAE 


Figure 3 
Also, it is not difficult to see that the method of solution presented 


for n = 2 earlier becomes increasingly complex as n increases, The 


2/ 


true computed mean output rate that is obtained by applying the 


approximate procedure, is found, as before, by determining the 


? ale 
Seed 
er 


combination of Ri +Rot ‘eine +R. =R that yields maximum R. However, 
it would seem that an iterative procedure could be developed to solve 
the problem efficiently for n = 3 or n = 4 on todays high speed 
computers. Efficient extension beyond this point would probably 
require knowledge the R,'s explicitly in terms of the service rates, 


queue lengths and flow distributions. 
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Pp, (Distribution from (3,1) to (4,1)) 


Ay» 42» M2143) = 0.5 A> faz = 0.4 
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PARAMETRIC STUDY 2. 
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P} (Distribution from (3,1) to (4,1)) 
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PARAMETRIC STUDY 3. 
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Abas and 4. (Second Stage Service Rates) 
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M34 and Ma5 (Third Stage Queue Length) 
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